Phantom phase power-law solution in f(G) gravity 
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Abstract: Power-law solutions for f(G) gravity coupled with perfect fluid have 
been studied for spatially flat universe. It is shown that despite the matter dominated 
and accelerating power-law solutions, the power-law solution exists for an special 
form of f{G) when this universe enters a Phantom phase. 
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I. INTRODUCTION 

Nowadays it is strongly believed that the universe is experiencing an accelerated ex- 
pansion. Recent observations from type la supernovae jl|] in associated with Large Scale 
Structure j^j] and Cosmic Microwave Background anisotropics jsj] have provided main ev- 
idence for this cosmic acceleration. These observations also suggest that our universe is 
spatially flat, and consists of about 70% dark energy (DE) with negative pressure, 30% dust 
matter (cold dark matter plus baryons), and negligible radiation. On the other hand the 
nature of dark energy is ambiguous. The simplest candidate of dark energy is a cosmolog- 
ical constant with the equation of state parameter w = —1. However, this scenario suffers 

n 

from serious problems like a huge fine tuning and the coincidence problem [4j. Alternative 
models of dark energy suggest a dynamical form of dark energy, which is often realized by 
one or two scalar fields. In this respect, dark energy has many dynamical components such 
as quintessence K-essence jg|, tachyon phantom ghost condensate and quintom 
9], and so forth. 

It is known that Einsteins theory of gravity may not describe gravity at very high ener- 



gies. The simplest alternative to general relativity is Brans-Dicke scalar-tensor theory 10|. 
Modified gravity also provides the natural gravitational alternative for dark energy 11]. 
Moreover, thanks to the different roles of gravitational terms relevant at small and at large 
curvature, the modified gravity presents natural unification of the early-time inflation and 
late-time acceleration . It may naturally describe the transition from non-phantom phase 
to phantom one without necessity to introduce the exotic matter. But among the most 
popular modified gravities which may successfully describe the cosmic speed-up is f(R) 



gravity. Very simple versions of such theory like 4 12] and 4 + R 2 13] may lead to the ef- 



fective quintessence/phantom late-time universe (to see solar system constraints on modified 



dark energy models refer to 



a 



). Another theory proposed as gravitational dark energy is 
scalar- Gauss- Bonnet gravity [15J which is closely related with the low-energy string effective 
action. In this proposal, the current acceleration of the universe may be caused by mixture 
of scalar phantom and (or) potential/stringy effects. The coexistence of matter dominated 
and accelerating power law solutions for this theory has already been shown In this 

paper, we extend these results to show the existence of Phantom phase power law solutions 
for an special form of f(G) gravity. 
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II. FIELD EQUATIONS FOR [R + f(G) + L m \ GRAVITY 



As an alternative to the f(R) action considered in Ref.|17j we consider the following f(G) 
action which describes Einstein's gravity coupled with perfect fluid plus a function of the 

n n 

Gauss-Bonnet term [181 ] . |19| 

1 



(1) 



S — J d^Xy/^g 

where k 2 = 8ttGn and the Gauss-Bonnet invariant is defined as follows 

G = R 2 - AR^Ri" + R^ Xa R^ Xa . (2) 
The field equations are obtained by varying the action with respect to g pu 

= ^(-R" U + l^ u R) + T>» + \g^f{G) - VgRR» v + 4f G R$R vp 

-2f G R^ T R; ar - AfaRWRpr + 2{V^ f G )R - 2g^(V 2 f G )R - 4(V p Wf G )R» p 

- 4(V p V7g)F + 4(V 2 f G )R^ + ig^(V p V a f G )R^ - 4(V P VJ G )R^ , (3) 

where f G = f'(G) and f GG = f"(G). The usual spatially-flat metric of Friedmann- 
Robert son- Walker (FRW) universe is chosen in agreement with observations 

3 

ds 2 = -dt 2 + a(t) 2 ^(rfx 4 ) 2 , (4) 
i=i 

where a(t) is the scale factor as a one-parameter function of the cosmological time t. Using 
this metric in the field equations fl3]) one obtains the first FRW equation 

- - 2 E 2 + Gfo - f(G) - 2AGH 3 f GG + Pm = , (5) 



where . denotes derivative with respect to time t and Hubble parameter H is defined by 
H = a/a. In the FRW universe, the energy conservation law can be expressed as the 
standard continuity equation 

p m + 3H (p m +p) = p m + 3H(1 + w)p rn = , (6) 

where p m is the matter energy density and p = wp m is the equation of state relating pressure 
p with energy density. From the continuity equation we obtain 

Pm(t) = P0t~ 3m(1+W) ■ (7) 



i 



For the metric (J4]), the Gauss-Bonnet invariant G and the Ricci scalar R may be defined as 
functions of the Hubble parameter 



G = 24(HH 2 + if 4 ), R = 6(H + 2H 2 



(8) 



III. EXACT MATTER DOMINANT AND ACCELERATING POWER-LAW 

SOLUTIONS 

We now assume an exact power-law solution for the field equations 



a(t) = a t m , 



(9) 



where m is a positive real number. Using this assumption in ([8]) leads us to the following 
results 



G 



24 

¥ 



m 3 (m — f 



G = -^m 3 (m - 1), 

6m . , 
R=— (2m -1). 

By substituting ([7]), f lTUj) and f lTTj) into ()5]) we obtain the Friedmann equation 



-G 2 }'gg + G/g — /g 



Gfi/2 



3m 



1/2 



Po 



(10) 

(11) 

(12) 
0. (13) 



m — 1 k 2 \8(m — l)y \24m 3 (m — 1) 

This is a differential equation for the function f(G) in G space. The general solution of this 
equation is obtained as 



/ (G) = CxG + C 2 G~^ m - 1) - 



1 



6m(m — 1 



G 1 ? + A mw G-* mil+w) 



where 



A, 



k A (m + l) 2 
po (m - 1) (13824m 9 (m - i) 3 )-> (1+u,) 



(14) 



(15) 



4 + m [3m(w + l)(w + 4/3) - Ihw - 19] ' 
and Ci, C2 are arbitrary constants of integration. This solution is in agreement with the one 



obtained in 



16[ and, as is explained there, we can without any loss of generality assume the 



constants C\ — Ci — 0. Hence, the required form of the function f(G) becomes 

f(G): 



1 

2 



6m(m — 1) 



k A {m + 1) 



+ A mw Gr m ^ 1+w ) 



(16) 
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First, we note that a real valued solution for f(G) requires the values m < or m > 1. 
While the former leads to a contracting universe and also causes a divergence at m — — 1 via 
the first term in the bracket, the latter is cosmo logically desirable for an expanding universe 
as follows. 

The case m — 1 leads to G = and i? = 6t~ 2 which is the general relativity limit with 
the power-law solution 

a(t) = a t, (17) 

and the energy density 

Pm(t)=P0t- 3(1+W) . (18) 

According to (J7J), it is easy to see that m = 2/[3(l + w)} indicates the reduction to general 
relativity, so for m = 1 the equation of state parameter is fixed by w = —1/3 which accounts 
for a negative pressure but not yet an accelerating universe. 

The case m > 1, leads to a nonzero real Gauss-Bonnet term G and a positive Ricci scalar 
R. However, in order to avoid divergence in the Gauss-Bonnet term we have to keep m and 
w away from the values for which A mw diverges according to the following equation 

4 + m [3m{w + l)(w + 4/3) - I5w - 19] = 0. (19) 

This case with m > 1 predicts an accelerating universe. Thus, power-law solutions of the 
type a(t) = a t m or H = ™ exist for the actions of the type [R + f(G) + L m ] with f(G) 
given by ( 1T41) except for those values of m which satisfy ( fT9i) . 



IV. EXACT PHANTOM PHASE POWER-LAW SOLUTION 

One may also study the power-law solutions where the universe enters a phantom phase 
leading to a Big Rip singularity. For this case, the general class of Hubble parameters and 
cosmological solutions are defined as 

HI 

H(t) = — f (20) 

a(t)=a (t s -t)- m , (21) 

where t s is the so called "Rip time" at future singularity. Again, using the above solution 
and repeating the similar calculations we obtain the following results 

p m (t) = Po t 3m{l+w \ (22) 
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24 

G = - — x rrf{m + 1), 



G 



R 



96 
{t. - 1) 5 
Qm 



m 3 (m + 1), 
(2m + 1). 



(t s -ty 

Substituting ([22]), fl23) and (J21J into the first FRW equation (ED we obtain 



(23) 
(24) 
(25) 



4 r 2 f 4-rr f ° 1/2 1 

■(j JGG + ^JG — JG- 



1/2 



Po 



G 



-§m(l+u/) 



0. (26) 



m+1 ■ " w 1 K 2 ^8(m + l)y r V 24m3 ( m + 1 ). 

This equation is easily recovered by the map m — > —m in the previous equation (fT3"|) . 
Therefore its solution is obtained by using the same map in (I14p as 



f(G) = C 1 G + C 2 G^ m+1) - - 



A; 4 (m-1) 2 + mu> 



where 



-A? 



8p (m + 1) (13824m 9 (m + l) 3 ) * m(1+ " } 



(27) 



(28) 



4 + m [3m(u; + l)(w + 4/3) + 15w + 19] ' 
Similar to the solutions in the previous section, we assume C\ — C' 2 — 0. Then, the required 
form of the function f{G) becomes 



/(G) = - 5 



'6m(m + 1) 
k\m - l) 2 



G^ + A n-z m ( 1 + w ) 



(29) 



Actually, m > leads to a real valued function f(G) according to ([2"7]) . However, demanding 
a Big Rip during the phantom phase, as the cosmic time t approaches t s , requires m > 1 
in ( 12T|) . However, m = 1 causes a divergence in /(G) through the first term of (|27|) in the 
bracket. Moreover, the Gauss-Bonnet term diverges through A mw for those values of m for 
which the following equation is satisfied 



4 + m [3m(w + l)(w + 4/3) + 15w + 19] = 0. 



(30) 



Therefore, power-law solutions in the phantom phase of the type a(t) = ao(t s — t)~ m or 
H{t) = Yzr t exist for the actions of the type [R + /(G) +L m ] with /(G) given by ( 1271) except 
for those values of m which satisfy fl30|) . 
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v. THE STABILITY ISSUE 



The stability issue of a large class of modified gravitational models has been discussed 
with particular emphasis to de Sitter solutions 20 2s| , [^ol , l^^] l^sl , 29 1. In the present 



modified gravity, namely [R+ f(G)], the stability issue leads to the following conditions [22 1 



Gofo — fo = 6if ( 



o- 



1 < 



d3 f" ' 



(31) 
(32) 



where the critical points are defined by 



R = 12Hq, Gq = 2AHq. 



(33) 



and /o, /o = (df /dG) , f = (d 2 f/dG 2 ) , H are suitable constants corresponding to the de 
Sitter solutions. In the case of Phantom phase power law solution, namely ( |29|) . the first 
condition reads as 



1 / 6m(m + 1) 
2~\i k\m-l) 2 



x 



1 - -m(l + u) ) A mw 



G | m(1+U)) _i = ^ 



(34) 



which implies that 



1 - -m(l+w) ) A mw > 0. 



By using (1351) . the second condition (1321 reads as 



, ,o /9 / 6m(m + 1) 



■-mil + uj'' 
4 



-m(l + u) + 1 J f 1 - -m(l + w) 



x 



\/6 



1 / 6m(m + r 

2 V k\m- iy 



< 9. 



(35) 



(36) 



Then, the model is stable around de Sitter solution if the arbitrary parameters also satisfy 
both the conditions (|35|) and (|36|) . 



s 



VI. CONCLUSION 

In the present paper we have considered an f(G) action which describes Einstein's gravity 
plus a function of the Gauss-Bonnet term. Then, by considering an exact power-law solution 
for the field equations we have obtained the Friedmann equation in spatially flat universe. 
The Friedmann equation appears as a differential equation for function f(G). We could 
obtain the solution of this equation and show that our model with this solution for f(G) 
has power-law solution of the type a(t) = a^t 171 expect for those values of m for which f{G) 



diverges. These solutions are in agreement with those obtained in [16[ . We have also studied 
the power-law solutions when the universe enters a Phantom phase. By considering such 
power-law solution for the field equations we have obtained the corresponding Friedmann 
equation. The solution f(G) of this differential equation is obtained and it is shown that 
the power-law solution in the phantom phase of the type a(t) = ao(t s — t)~ m exists for this 
f(G) except for those values of m for which this function diverges. 
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